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We introduce two models of quasi-discrete optical systems: an array of waveguides doped by 
four-level N-type atoms, and a nonlinear checkerboard pattern, formed by doping with three-level 
atoms of the A- type. The dopant atoms are driven by external fields, to induce the effect of the 
electromagnetically-induced transparency (EIT). These active systems offer advantages and addition 
degrees of freedom, in comparison with ordinary passive waveguiding systems. In the array of active 
waveguides, the driving field may adjust linear and nonlinear propagation regimes for a probe signal. 
The nonlinear checkerboard system supports the transmission of stable spatial solitons and their 
"fuzzy" counterparts, straight or oblique. 

PACS numbers: 42.65.Tg, 42.82.Et, 42.50.Gy, 42.50.Ct 



I. INTRODUCTION 

The transmission of light in discrete arrays of evanescent ly coupled waveguides is a topic of great interest in optics. 
The arrays are prime examples of systems in which the discrete optical dynamics can be observed and investigated 
Optical fields propagating in such settings exhibit a great number of novel phenomena ^^\. However, traditional 
coupled arrays or lattices, such as arrays of waveguides made of AlGaAs |^ or periodically poled lithium niobate 
(PPLN) [5], virtual lattices created in photorefractive crystals (PhCs) [g^ and liquid crystals [7] by means of the 
optical induction, etc., are built of passive ingredients. 

On the other hand, it is well known that active elements, such as atoms with a near-resonant transition frequency, 
may lend the medium a number of specific optical characteristics, such as strong dispersion, a complex dielectric 
constant, and the strong variation of the dispersion relation near the resonance. Therefore, waveguide arrays made 
of active elements may offer low thresholds, in comparison with their passive counterparts, and possibilities for the 
"management" of their waveguiding characteristics. In the ID (one-dimensional) case, active discrete systems were 
previously studied in detail in the form of resonantly-absorbing Bragg reflectors (RABRs) [8], which are used to 
demonstrate optical switching [9], storage [10], and nonlinear conversion 

Recently, a 2D (two-dimensional) "imaginary-part photonic crystal" (IPPhC, i.e., a medium with a periodic varia- 
tion of the imaginary part of the refractive index) is realized by means of the techniques of multi-beam-interference 
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holography, hthography and back-fihing [l2^, which ahows one to create a spatiahy-structured distribution of the 
active material. For example, active substance Rhodamine B can be doped into the homogeneous SU8 background 
to form an IPPhC. In this structure, the real part of the refractive index is constant if the probe wavelength is far 
detuned from the resonance. However, the imaginary part of refractive index affects the real part, which becomes 
a conspicuous effect close to the resonance. Thus, in the vicinity the absorption window, the IPPhC also acts as a 
traditional PhC. Very recently, a laser system in a medium featuring a periodic distribution of loss, which is akin to 
IPPhC, is demonstrated in the experiment [isl ]. 

In this work we propose two new varieties of active light-guiding systems. In Section II, we show that the introduction 
of an active material into the PhC provides for a way to create active structures in the form of a coupled waveguide 
arrays. The difference of this system from the RABR is that guides light not across the periodic structure, but rather 
along it. In Section III, we introduce a checkerboard system, which is built of alternating linear and the nonlinear 
square cells in the x-z plane. These systems may be controlled ("managed") via the effect of the electromagnetic- 
induced transparency (EIT). 



II. WAVEGUIDING ARRAYS CONTROLLED BY THE ELECTROMAGNETICALLY-INDUCED 

TRANSPARENCY 

In this section we consider the possibility to use N-type near-resonant four- level atoms as the active dopant. The 
respective scheme of the energy levels is shown in FIG. 1(a), where |1) and |2) are the ground and a metastable states, 
respectively. These two states have the same parity of their wave functions, which is opposite to that of states |3) and 
|4>. 

As a part of the scheme, we assume that a weak probe wave, Ep^ with Rabi frequency ftp = psiEp/h is acting 
on transition |1) |3), with single-photon detuning Ai. Here psi (which is assumed real) is the matrix element of 
the dipole transition between |1) and |3) . Further, a traveling- wave field with Rabi frequency Qc drives the atomic 
transition |2) |3) with detuning Ac = Ai, hence the two-photon detuning is given by (5 = Ai — Ac = 0. As 
another ingredient of the EIT scheme, an optical-induction field with Rabi frequency Qs induces transition |2) |4), 
with detuning A2. The decay rate for level |n) is 7^. Here we neglect 71 and 72, and assume 73 ?^ 74 = 7. 

The Hamiltonian of the system is: 

H =^h^i\l){l\--[Qpe-^--'\3){l\^ 
1=1 

l^ce"'^'''|3)(2| + nse-'^''\^){2\ + H.C]. (1) 

where uji is the eigenfrequency of the l-th level. The EIT effect means that, when the probe is exactly at the two- 
photon resonance (5 = 0), and the atoms are prepared in the ground state, the linear absorption of the probe vanishes, 
irrespective of the single-photon detuning [14]. 

The quasi-discrete system, which is introduced in this section, is illustrated by panel (bl) in FIG. 1, which shows 
the distribution of the concentration of the active component in the ID case. The transverse width of the waveguide 
is di, the interval between the waveguides is ^2, and Nq is the density of the N-type atoms inside the waveguides. 
Coefficients are chosen as per experimental in the case of Y2Si05 doped with Pr^+ (Pr:YSO) [15]: the density of 
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FIG. 1: (Color online) (a) The energy- level diagram of the N-type atom, (bl) The structure of the active waveguide array: the 
transverse width of the waveguide is di , the interval between the waveguides is d2, and the density of the active atoms in the 
waveguide is A^o • (b2) The respective effective periodic potential with Vb = —Nopli\Qsf/'2eohA2\^cf • (c) The energy-level 
diagram of the A-type atom, (d) The checkerboard system: black and white square cells depict areas which are, respectively, 
doped with the active material, or are left undoped. 




active atoms inside the waveguides is A^o = 1-0 x lO^^cm^ (which corresponds to the dopant concentration ^0.1%), 
p3i = 1.18 X 10~^^C-m, and 7 = 30 kHz, the probe wavelength being 605 nm. Linear and the nonlinear properties of 
this system are detailed below. 

(1) Linear properties of the system 

If we set Ai = and assume A2 ^ ^c,7 , then the absorption of the probe can be neglected. One can easily find 
the one-step steady-state solution for the density-matrix element of the transition between |1) and |3) , cf. Ref. jicj : 

Therefore, the contribution of the resonant atoms to the polarization experienced by the probe is ^ = 2Npsipsi, 
cf. Ref. [17] (recall N is the dopant density). The (1+1)D paraxial propagation equation for the slowly varying 
envelope of the probe field, Ep, is 

2zkp^Ep = -^Ep-^^, (3) 

where kp = 27rn/Ap is the wavenumber of the probe, and n the background refractive index. The substitution of Eq. 
(2) into the last term of Eq. (3) yields a scaled linear Schrodinger equation, 

4^ = -^^f/+n.)t/, (4) 

where ( = kpZ, r] = kpX, U = flp/j, and the effective potential, 

is induced by the concentration distribution N{r]) , as shown in panel (bl) of FIG. 1. This potential is induced by 
the giant Kerr effect controlled by the Rabi frequency Qs [iSl- In fact, this potential emulates a difference in the 
local refractive index contrast between the active (doped) and passive (undoped) regions. For A2 > (here, we 
choose A2 = IOO7), the shape the potential is shown in panel (b2) of FIG. 1, which is similar to the Kronig-Penney 
potentials corresponding to the tight-binding model in solid-state physics . The depth of local wells in the periodic 
potential is defined by the intensity ratio, l^l^p/l^lcP, the corresponding refractive-index contrast between the active 
and passive regions being An = p^iNo\Qs\'^ /'^^o^^2\^c\'^ 
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FIG. 2: (Color online) Simulations of the linear transmission of light in the arrays of active waveguides. Here, we chose 
di = d2 = 8, Qc = 27 and A2 = IOO7. The incident beam (at ( = 0) is [/ = Aq e^p{-r]^ /df) with Aq = 0.07. (a) 
\Qsf = 0.036|Qc|^ (i.e., An = 0.09). (b) \Qsf = 0.064\Qcf (i.e.. An = 0.16). (c) l^^^l^ = 0.100|Qc|^ (i.e.. An = 0.25). (d) 
\ns\^ = 0.144|Qc|^ (i.e.. An = 0.36). 



Figure 2 displays the propagation of the probe in such a system at different values of l^l^p/ir^cP (i-e. different 
values of An), produced by numerical simulations of Eq. (4). This figure demonstrates that the quasi-discrete 
diffraction naturally gets suppressed when the depth of the potential, i.e., I^l^p/ir^cp or An, increases, resulting in a 
reduced coupling between local waveguides. Thus, the diffraction in the present setting may be efficiently controlled 
by varying the magnitude of l^l^p/l^lcP- 

(2) Nonlinear properties of the system 

If both detunings satisfy conditions Ai, A2 ^ ^c, 7, then the absorption of the probe (the losses) may be neglected. 
Moreover Ai 7^ produces an enhanced Kerr nonlinearity [20], provided that the density matrix element pai, which 
accounts for transitions between |1) and |3), is taken into regard, to describe the third-order effect. The steady-state 
solution for this matrix element can be obtained in the following from, cf. Ref. (2l|: 



(1) , (2) , (3) _ o l^^l^ o ra\ 

With regard to this result and the above definition, U = Eq. (3) changes its form into that of the standard 
nonlinear Schrodinger (NLS) equation [22]: 

= " ^ ^^''^^ ^ (7) 
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where the effective nonhnear coefficient is 

Thus, is affected by the modifications of the refractive index of two types: (1) a periodic change of the linear 
index induced by Qs via the giant Kerr effect; (2) the nonlinear change under the action of Ep itself, via the enhanced 
self-Kerr effect. The sign of detuning Ai determines whether the latter effect gives rise to the self-focusing (Ai < 0) 
or self-defocusing (Ai > 0) sign of the nonlinearity. When the tunneling coupling between adjacent waveguides is 
balanced by the nonlinearity, quasi-discrete solitons 12311 can be formed in this system, which is similar to those in the 
traditional coupled arrays made of passive materials [2^. 

III. SOLITONS IN THE CHECKERBOARD SYSTEM CONTROLLED BY THE 
ELECTROMAGNETICALLY-INDUCED TRANSPARENCY 

In this section, we assume that the active material is doped periodically in both the transverse and propagation 
directions (i.e., along the x and z axes, respectively). The corresponding density distribution of the active material 
is N{x,z) = NoR{x,z), where R{x^z) is a dimensionless structural function of the distribution. Here, we adopt for 
R{x, z) the form checkerboard form depicted in FIG. 1(d). The white cells, with i?(x, z) = , are areas which are not 
subject to the doping, while black cells, with R{x^ z) = 1 , depict areas doped by the active material. 

Note that the formation of 2D spatial solitons in the checkerboard-shaped linear potential was considered in Ref. 
[25^ . assuming that the probe beam was shone along the uniform direction in the bulk medium equipped with the 
transverse checkerboard structure. Here, the difference is that the modulation is applied to the nonlinear term, and 
light propagated across the structure. 

If we turn off the control field Vts^ the four- level N-type atomic configuration reduces to the three- level one of 
the A- type, see FIG. 1(c). In this case, the linear-refractive-index contrast between the active and the passive areas 
vanishes, which leaves only the nonlinearity modulation in action. 

Again assuming that the detuning is much larger than the decay rate, Ai 7 , the absorption of the probe may 
be ignored as before. Accordingly, in the present setting Eq. (6) is rewritten as: 

...^.S'+.g'+.S'-^Op (9) 

and Eq. (3) changes into 

i2 



where we define 



Therefore, the white and black cells, with hz{r]X) = and n{ri^Q = act as linear and nonlinear elements, 
respectively. 

The periodic modulation of the nonlinearity in Eq. (10) places this equation into the broad class of models 
with nonlinear lattices (see original works [26^ and review [27]). However, the present checkerboard pattern of the 
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modulation in the longitudinal and transverse direction was not studied in previous works. More general patterns, that 
should be studied separately, may be represented by arrays of isolated black squares set against the white background, 
or vice versa. 

Below, we choose values Ai = — IOO7, = 0-57, with the other parameters taken as before. This yields the 
nonhnearity- modulation amplitude hzo = —10 . The probe field is launched at C = as a Gaussian 

U{7i,0) = Aexp[-{7j-af/W% (12) 

with the central point at = a. For the simulations, we take the checkerboard with squares of size 8 x 8, as shown 
in FIG. 1(d). 

The simulations of the evolution of the Gaussian in the framework of Eq. (10) were carried out by means of the 
split-step Fourier method. First, we chose the amplitude and width of Gaussian (12) as, ^4=0. 075, 0.065, 0.055 and 
W=S^ with the center placed at the mid-point of the nonlinear cell (a = 4). Results of the simulations are displayed 
in FIG. 3(a)-(c). In particular, FIG. 3(a) shows that the probe field with A=0.075 propagates without decay and 
distortion over the distance longer than z = 100 x W^^ i.e., ~100 diffraction lengths. This dynamical regime may be 
naturally identified as a stable soliton. On the other hand, in FIG. 3 (b), with the Gaussian's amplitude ^4=0. 065, 
the probe field forms a fuzzy beam, which, nevertheless, avoids decay over the distance exceeding z = 100 x W^. 

However, for ^4=0. 055, FIG. 3(c) demonstrates that the input cannot form a robust beam and rapidly decays. 
Therefore, there must be internal borders (i.e. thresholds) separating the stable solitons, fuzzy beams, and decaying 
ones, in the plane of the width and amplitude of the Gaussian inputs. These borders are plotted in FIG. 3(d). 

Further simulations, displayed in FIG. 4(a), show that stable straight solitons can be formed as well if the Gaussian 
is launched at the midpoint of the linear cell. However, if the center of the Gaussian does not coincide with the center 
of the linear or nonlinear cell, the propagation of the soliton beam becomes oblique, see an example in FIG. 4(b). 

It may be interesting to consider the oblique propagation of beams across the checkerboard, induced by the ap- 
plication of a lateral kick to the input. Another issue of obvious interest is the interaction of beams in this setting. 
These generalizations will be reported elsewhere. 



IV. CONCLUSION 



In this work, we have proposed a method for building systems of coupled active waveguides, and also a new system 
with the checkerboard pattern of the modulation of the nonhnearity coefficient. These settings can be created using the 
appropriate doping patterns of the N-type four- level and A- type three- level resonant atoms, respectively, and driving 
them by means of the EIT mechanism. Such active systems offer certain advantages, admitting more possibilities for 
the design and management, in comparison to the passive media. Firstly, the nonhnearity in the active systems can 
be switched via the sign of detuning Ai: the incident beam with Ai < and Ai > will experience the action of 
the self- focusing and self-defocusing nonhnearity, respectively. Next, it is well known that the EIT may be efficiently 
applied to few-photon settings, especially in the nonlinear regimen Q, Therefore, the systems introduced here 
may, in principle, offer an advantage for handling quantum and non-classical light beams, composed of few photons. 
Further, it is well known that the group velocity of the probe can be coherent controlled j29| and tuned to a very small 
value under the action of the EIT ^] , which implies that the probe signal in this system can be trapped in the from 
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FIG. 3: (Color online) Simulations of the evolution of the Gaussian input probe with width = 8 in the checkerboard model, 
(a) With amplitude A=0.075, the probe propagates in a stable fashion over the distance of L = 200 x = 12800. This beam 
may be called a stable soliton. (b) With A = 0.065, the probe propagates, keeping a fuzzy shape, still featuring robustness 
against the diffractive spreading-out. (c) With A = 0.055, the probe quickly spreads out. (d) In the plane of the width and 
amplitude if the input Gaussian, A^iW) is the border between the stable and fuzzy beams, while A^iW) is the border between 
the fuzzy but robust beams and decaying inputs. 
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FIG. 4: (Color online) (a) The probe with parameters A = 0.075, = 8, and a = —4 is launched at the mid-point of the 
linear cell. The probe propagates in a stable fashion over distance L = 200 x = 12800, therefore it is categorized as a stable 
soliton. (b) The probe with A = 0.075, W — 8, a = propagate obliquely if it is launched off the mid-point of a linear or 
nonlinear cell. 



of the slow light. Thus, various applications of the EIT, such as dark-state polar itons [sil, the few-photon four- wave 
mixing ultra- weak and ultraslow light [sS], etc., may be realized in the systems proposed here. Furthermore, 
using properly designed holographic patterns, various complex spatial structures of the distribution of the dopant 
concentration can be photoinduced in the 2D geometry, such as quasi-crystals [34], honeycomb lattices [35], defect 



lattices [36], ring lattices 



etc., in addition to the simplest checkerboard patterns analyzed herein. Such 2D 



structures may have their own spectrum of potential applications. 
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